8. REFRACCION SISMICA 


8.1 DESCRIPCION GENERAL 

Dentro de los metodos sismicos de la geofisica aplicada se encuentran los de 
refraccion y reflexion sismica. En estos metodos se mide el tiempo de 
propagacion de las ondas elasticas, transcurrido entre un sitio donde se generan 
ondas sismicas y la llegada de estas a diferentes puntos de observacion. Para 
esto se disponen una serie de sensores en linea recta a distancias conocidas, 
formando lo que se conoce como tendido slsmico o Ifnea de refraccion - o 
reflexion - sismica. 

A una distancia conocida del extremo del tendido, en el punto de disparo, se 
generan ondas sismicas, - con la ayuda de un martillo o por la detonacion de 
explosivos las cuales inducen vibraciones en el terreno que son detectadas por 
cada uno de los sensores en el tendido. 

El equipo basico consiste de los sensores; la unidad de adquisicion, en donde 
se almacenan los movimientos del terreno detectados por cada sensor; los 
cables de conexion entre los sensores y la unidad de adquisicion; el cable del 
trigger, que se encarga de marcar el momento de inicio de registro en la unidad 
de adquisicion. 
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Los registros de cada sensor tienen informacion de los movimientos del terreno 
en funcion del tiempo y son conocidos como sismogramas. Estos son 
analizados en la refraccion sismica para obtener el tiempo de llegada de las 
primeras ondas a cada sensor desde el punto de disparo, y en la reflexion para 
obtener informacion de las ondas que son reflejadas en las diferentes interfaces 
de suelo, para lo cual es estudiado el sismograma completo. 

Como se vera mas adelante una de las aplicaciones del metodo de refraccion 
sismica en la ingenieria civil es la determinacion de la profundidad al basamento. 
Con este fin fue aplicado este metodo en este trabajo de grado en la zona de 
tesis. 

8.1.1 Aplicaciones en la Ingenieria Civil. 

La aplicacion mas comun de la refraccion sismica en la ingenieria civil es para la 
determinacion de la profundidad a basamento en los proyectos de construccion 
de represas y grandes hidroelectricas, y para la determinacion de las 
condiciones (meteorizacion, fracturacion) y competencia de la roca en donde se 
asentaran las estructuras, asi como por donde se realizaran los tuneles. 
Tambien es muy util para deteccion de fallas geologicas. 

En el caso de contextos urbanos la refraccion resulta util para la determinacion 
de la profundidad a basamento y el perfil de velocidades de onda P y S; y para la 
extrapolacion lateral de perforaciones puntuales de suelos. 

El metodo utiliza la llegada de las primeras ondas a los geofonos, ondas P, pero 
como tambien las llegadas de las ondas S, de tal manera que se pueden 
determinar la relacion de Poisson y otros modulos dinamicos. 

Sarria (1996) enuncia otras potenciales aplicaciones del metodo: a) utilizacion 
del ruido sismico para determinar el modulo dinamico G en masas de arcilla; b) 
evaluacion del amortiguamiento; c) evaluacion de los limites de Atterberg; d) 
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determinacion de modulos E y G en mecanica de rocas; e) explotacion de 
canteras; f) ubicacion de sondeos en roca y g) para determinar la capacidad de 
carga de los pilotes. 

8.1.2 Alcances y limitaciones del metodo. 

En terminos de la Ingenieria Civil, y el estudio dinamico de los suelos los 
alcances y limitaciones del metodo serian los siguientes: 

Alcances. 

• Detecta variaciones tanto en profundidad como en la horizontal de la 
velocidad de la onda P ( y de la S). 

• Permite la deteccion de la profundidad a basamento y de su relieve, 
dependiendo de variables como longitud del tendido, energia de la fuente 
sismica, velocidades de los suelos. 

Limitaciones. 

• Solo funciona cuando la velocidad de propagacion de las ondas aumenta con 
la profundidad. En el caso de suelos con capas intermedias de menor 
velocidad el metodo arrojaria resultados erroneos. 

• Para el caso de aplicaciones urbanas de la Ingenieria Civil, el Metodo de 
Refraccion Sismica esta limitado por la disponibilidad de zonas descubiertas 
con suficiente extension. La longitud del tendido en superficie esta 
directamente relacionada con el alcance de la exploracion en profundidad. 

8.2 PROPAGACION Y TRAYECTORIA DE LAS ONDAS 

Cuando se generan ondas sismicas, a partir de golpes en el suelo con una porra, 
o con explosiones de polvora, estas incluyen tanto ondas sismicas internas, - 
Primarias y Secundarias -, como superficiales ondas - Love y Rayleigh -. 
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Las ondas P, tambien conocidas como ondas longitudinales, son las de mayor 
interes en la refraccion sismica. 

Las leyes que rigen la propagacion y la trayectoria de las ondas sismicas en la 
refraccion, son las mismas que se utilizan en optica: 

(1) Principio de Huygens. 

(2) Principio de Fermat, y 

(3) Ley de refraccion (o de Snell), 

las cuales fueron explicadas en el Marco Teorico, Capitulo 3. A continuacion 
solamente se retoma y amplia el concepto de la Ley de refraccion. 

Ley de refraccion. 

Como consecuencia del Principio de Huygens y/o del principio de Fermat, la Ley 
de refraccion dice que el seno del angulo incidente es al seno del angulo de 
refraccion como la velocidad de la onda incidente es a la velocidad de la 
correspondiente onda refractada. 

Para explicar la trayectoria de las ondas en el metodo de la Refraccion sismica, 
consideremos un medio, con velocidad Cl, que suprayace un medio semi- 
infinito, con velocidad C2, mayor que Cl (Figura 31). Una vez se han generado 
las ondas en el punto de disparo, estas empiezan a viajar por el medio superior 
conformando unos frentes de onda en el espacio. 
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Figura 31 . Modelo de dos capas, la inferior de mayor velocidad 

Al hacer un corte vertical por el punto de disparo, el frente de ondas luciria como 
se ilustra (Figura 32-a). Dicho frente se conocen como frente de ondas directas. 
En la parte b) de la Figura 32 el frente de ondas se ha encontrado con el llmite 
de los medios y ocurren las primeras refracciones hacia la capa inferior. En la 
parte c), ha pasado mas tiempo y se pueden observar claramente 3 frentes de 
onda: 1. de las ondas directas; 2. de las ondas refractadas hacia la capa inferior, 
y 3. de las reflejadas hacia la capa superior. Al observar en detalle puede 
identificarse un cuarto frente de ondas. El frente de ondas refractado hacia la 
capa inferior, no tiene una curvatura constante, de tal manera que corresponde a 
dos frentes de onda, el que se refracta hacia abajo, y el que se refracta hacia la 
capa superior. Como se puede observar, este frente de ondas esta mas alejado 
del punto disparo que el frente de ondas directas en la primera capa, por lo que 
Negara mas rapido a los geofonos donde aun no habia llegado el frente de ondas 
directas. 

En la parte d), de la figura, ha pasado aun mas tiempo desde el momento de 
disparo, y los 4 frentes de onda se diferencian claramente. 
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Figura 32. Propagation de las ondas en un medio de dos capas 


El frente de ondas refractadas hacia el medio superior se genera cuando los 
rayos provenientes de la fuente alcanzan en angulo critico, i c , la interfaz entre 

los medios. Como se explico en el Marco Teorico, Capitulo 3, el angulo 
refractado tiene 90° con respecto a la normal, de tal manera que el 

sen(i 2 ) = sen{ 90) = 1 . 

La refraccion a 90° del angulo critico, ilustrada en la Figura 33, implica que las 
ondas no se propagan por la capa inferior, sino por el contacto entra ambas 
capas - es decir por la superficie de refraccion - con la velocidad v2 de la capa 
inferior, siguiendo la ley de los recorridos rmnimos o Principio de Fermat (p. ej. 
Boyd, 1999; Masuda, 1983). 
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Adaptado de Boyd, 1999. 
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Figura 33. Angulo critico de refraccion i c . 


Fenomenos en la propagacion. 

Cuando el medio en que se propagan las ondas sismicas no es homogeneo, se 
producen los fenomenos de difraccion, dispersion y scattering (p. ej. Cantos, 
1973). 

Difraccion. Desvio de los rayos, en cierta extension, ocurrido cuando se limita 
parte del frente de ondas (e. g. Tipler, 1990). 

Dispersion. Es la variacion de la velocidad de una onda con el cambio de 
frecuencia. En un medio elastico homogeneo no hay dispersion, pero si la hay en 
un medio imperfectamente elastico como en la tierra. En refraccion sismica no 
hay evidencia de que exista dispersion apreciable (p. ej. Cantos, 1973), excepto 
cuando se usan explosivos en inmediaciones de la explosion. 

Scattering. Corresponde a la formacion de pequenas ondas que propagan la 
energia en todas las direcciones. Se produce cuando un frente de ondas choca 
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con particulas libres u objetos pequenos comparados con su longitud de onda. 
Este fenomeno no es mayor para frecuencias altas. Parte de lo que se considera 
“ruido” en un registro puede deberse a este fenomeno ya que produce energia 
distribuida al azaren superficie (p. ej. Cantos, 1973). 


La disminucion de la energia sismica con la distancia, causada por los tres 
fenomenos explicados anteriormente, va acompanada de perdidas debidas a la 
absorcion de la energia, produciendo amortiguamiento. Cuando el impulso 
sismico viaja a traves de las diferentes capas las altas frecuencias son 
absorbidas mas rapidamente que las bajas frecuencias. 

8.3 EQUIPO DE MEDICION 

Las componentes del equipo de medicion de refraccion sismica son (Figura 34): 

• Fuente de generacion de ondas sismicas. 

• Deteccion de los movimientos del terreno. 

• Adquisicion y almacenamiento. 
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Figura 34. Esquema general del equipo de refraccion sismica. 


Fuente de generacion de ondas. 

Las ondas que se utilizan en refraccion son generadas por una perturbacion 
artificial instantanea, que se conoce como impulso sismico. Lo que se busca 
con este es generar el tipo de ondas sismicas (de volumen y superficiales), 
producidas por un unico “evento” de duracion “instantanea”, para que no haya 
superposicion de ondas (de diferentes eventos) en los movimientos del terreno 
detectados por los geofonos. 
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Para esto se buscan fuentes de generacion que se puedan controlar en 
terminos del tiempo de inicio y localization (p. ej. Jakosky, 1950; Cantos, 1980). 
Existen tres tipos: 

• Fuentes de impacto. Generalmente martillos o porras. Como la energia 
transmitida al suelo por este tipo de fuente no es muy grande, se apilan 
varias decenas de golpes para modelar mejor las llegadas y suprimir el ruido. 
Tambien se utilizan otros medios mecanicos, por ejemplo dejando caer un 
gran peso de una altura de 2 - 3 metros. La energia asociada con cada uno 
de los golpes depende la energia cinetica (E) que relaciona la masa del 
martillo (m) y la velocidad aplicada al martillo (v), obteniendo E = m*v 2 . 

• Cargas explosivas. De mayor energia, son usadas especialmente para 
prospeccion petrolera. La explosion puede ocurrir en un tiempo de micro a 
milisegundos, dependiendo de la naturaleza y cantidad del explosivo y, del 
material que rodea sitio de explosion (p. ej. Cantos, 1980). 

• Disparos. Se utilizan balas o cartuchos de fogueo. La energia es mayor que 
la generada por martillo. 

Deteccion de los movimientos del terreno. 

A traves de geofonos de una componente vertical, el movimiento del terreno es 
observado en diferentes puntos a lo largo del tendido de refraccion sismica. 
Actualmente se usan 12, 24 6 48 geofonos. Estos sensores exigen mayor 
resistencia mecanica 15 que aquellos usados en la sismologia tradicional (p. ej. 
Jakosky, 1950) debido a que en refraccion se requiere geofonos con frecuencias 
naturales de vibracion mucho mayores, entre 8 y 40 Hz. 

Adquisicion y almacenamiento. 

Las partes que componen el equipo de adquisicion y almacenamiento son: 


15 La resistencia mecanica depende de la masa y de la constante de rigidez del sistema 
mecanico del sensor. 
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• Unidad de apilamiento y digitalizacion. Corresponde a la unidad donde se 
adquiere, digitaliza y procesan los datos. Tiene puertos especializados para 
recibir las senales enviadas por los geofonos; a cada geofono le corresponde 
un canal y sus senales son filtradas y digitalizadas de acuerdo con las 
opciones definidas por el usuario. 

• Cables sfsmicos. Comunmente se usan dos cables para conectar en cada 
uno la mitad de los geofonos. Cada cable contiene a su vez cableado para 
llevar la serial de cada geofono a la unidad de adquisicion; cada uno de estos 
se llama canal. 

• “Trigger”. Cable conectado apropiadamente a la fuente sismica, de tal 
manera que en el instante en que se golpea el suelo con el martillo o cuando 
la carga explosiva es detonada, el sistema de registro empieza a grabar. 

8.4 CURVAS TIEMPO-DISTANCIA 

Las curvas tiempo distancia se construyen con los tiempos de llegada de las 
ondas P a cada uno de los sensores, y la distancia de cada sensor al punto de 
disparo. 

Los tiempos de llegada pueden ser leidos directamente en la pantalla de la 
unidad de adquisicion, (o de una impresion realizada desde dicha unidad); o en 
la computadora, a traves de un software especial izado, una vez los datos hayan 
sido transferidos de la unidad de adquisicion al computador. El conjunto de 
registros (del total de geofonos) como se pueden apreciar en una impresion 
realizada desde la unidad de adquisicion, son mostrados en la Figura 35. 

Para tener informacion mas detallada del subsuelo a analizar, se realizan 
tendidos de refraccion conjugados, llamados tendido directo, reverso e 
intermedio. En el primero el punto de disparo se ubica en un extremo del tendido 
a una distancia conocida, mientras que en el segundo el punto de disparo se 
ubica al otro extremo del tendido; en el tendido intermedio, el punto de disparo 
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es colocado usualmente hacia el centra del tendido. Cuando se van a estudiar 
terrenos muy extensos, o cuando se necesita informacion muy detallada, se 
acostumbra hacer tendidos traslapados que permiten modelar mejor la 
topografia de las discontinuidades. 
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Figura 35. Registro tipico de los sismogramas 


Las curvas de diferentes puntos de disparo (p. ej. directo, reverso) de un mismo 
tendido se dibujan sobre una misma grafica como se ilustra a continuacion 
(Figura 36). 



Figura 36. Curva tiempo - distancia para disparos directo y reverso 
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Cuando se realizan tendidos de refraccion traslapados, las curvas tiempo - 
distancia tambien se dibujan en relacion con el traslape de los tendidos. 

La distancia critica 16 x c , es aquella medida entre el punto disparo y el sitio donde 

emerge la primera onda refractada en superficie. Dependiendo de las 
velocidades de la capa superior, y del refractor - o capa inferior y de la 
profundidad a este, la distancia critica puede ser o no menor que la longitud del 
tendido en observacion. 

En la Figura 37 se muestra la relacion entre el contraste de velocidades de la 
primera capa vl y del refractor v2, versus la relacion entre la distancia critica x c y 
la profundidad A . 
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Figura 37. Relacion entre la distancia critica y la profundidad 


16 Conocida en ingles como cross over distance 
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A medida que aumenta la relacion v2/vl la relacion xjh disminuye. Como regia 

practica, las longitudes de los tendidos de refraccion deben ser mayores que el 
doble de la profundidad al refractor para observar refracciones sin interferencias 
indebidas de las ondas P originales (Sheriff & Geldart , 1991). 

8.5 SUPOSICIONES DEL METODO 

La parte de la fisica que trata los rayos como lineas es llamada optica 
geometrica, y esta basada en las siguientes suposiciones (p. ej. Masuda, 1981): 

a) Los rayos viajan en lineas rectas por el vacio o por un medio homogeneo. 

b) En el caso hipotetico que un rayo se devolviera, lo haria siguiendo la misma 
trayectoria hasta llegar al punto de origen. 

c) No hay interaccion entre rayos. Los rayos pueden avanzar 
independientemente. 

d) Los rayos siguen las leyes de la reflexion y refraccion. 

Las ondas sismicas se propagan como un conjunto de rayos a traves de las 
rocas, los suelos, o la litosfera en general. La suposicion general del metodo de 
refraccion sismica es: 

• Las velocidades aumentan con la profundidad: el metodo no permite 
identificar capas o estratos de suelo con velocidades inferiores a la superior 
(“inversiones de velocidad”). 

8.6 PRINCIPIOS GENERALES DE INTERPRET ACION EN REFRACCION 

Una vez construidas las curvas tiempo-distancia (t-x) se procede a su 
interpretacion. La tarea principal es identificar las secciones de las curvas que 
pertenecen a un mismo refractor. El conjunto de puntos que pertenecen a un 
refractor conforma lo que se conoce como dromocrona. La identificacion de las 
dromocronas es la parte mas importante de la interpretacion de los datos de 
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refraccion sismica. Hay algunas caracteristicas del subsuelo que se pueden 
prestar para malas interpretaciones, a saber: 

• Un cambio de pendiente del curva T-x no significa necesariamente un cambio 
de refractor, sino que puede significar un cambio de pendiente del primer 
refractor, (p. ej. Cantos, 1989). 

• Cuando existe un estrato o una capa delgada de suelo cuya velocidad es 
menor que la de la capa superior, no hay refraccion critica, de tal manera que 
no habria indicios de su presencia en las primeras llegadas en cada punto de 
la linea de sismica. ( p. ej. Sheriff & Geldart, 1 991 :282). 

• Cuando existe una capa demasiado delgada, a pesar de tener velocidades 
mayores no alcanza a producir primeros arribos por el hecho mismo de ser 
tan delgada (p. ej. Sheriff & Geldart, 1991:283, Sarria, 1996). 

Hay 5 principios (y/o leyes) generales que conforman la base para la 
interpretacion de un conjunto de datos de refraccion sismica. Dichos principios 
que se explican a continuacion: 

1. Leyes de Snell. 

2. Ley de las velocidades aparentes. 

3. Principios de reciprocidad. 

4. Principio del tiempo de intercepto en el origen. 

5. Principio de paralelismo. 

Debido a que el primero fue explicado ampliamente en el Marco Teorico Capitulo 
3 a continuacion se explican los otros cuatro. 

8.6.1 Ley de las velocidades aparentes. 

La ley de las velocidades aparentes dice que la velocidad con que aparenta 
transmitirse una onda en un cierto punto de la superficie del suelo es igual al 
cociente entre la velocidad superficial y el seno del angulo de emergencia, 
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tornados ambos en dicho punto. Donde el angulo de emergencia es formado por 
la onda emergente con la superficie (p. ej. Cantos, 1973). 


Sea un code vertical del terreno (Figura 38) y un frente de ondas GL llegando a 
la superficie en G . Considerese dos rayos infinitamente proximos llegando a dos 
geofonos G y G' separados una distancia Ax . Sean GA y G'B los tiempos de 
llegada. La velocidad con que aparente transmitirse los rayos G y G'sera 

V a = \ consideremos adicionalmente que el frente de ondas se desplaza de 

L a G’con su velocidad superficial V 0 , donde A^=F 0 y como As = Ax*sen(G 0 ) , 
se obtiene la velocidad aparente V a : 


V = ^ = 


As 


K 


con V a > V 0 


At sen (9 0 ) * At sen{ 9 0 ) 

La velocidad aparente en el punto A, sera la tangente a la dromocrona AB en el 

\x 

punto A, ya que V a = lim — , cuandoAt -> 0 teniendo finalmente que: 

At 


V„ = 


dx 

dt 


8.6.2 Principio de reciprocidad. 

Este principio establece que el tiempo de propagacion de una onda sismica del 
un punto A al B, es el mismo que el de B hacia A. Esto es una consecuencia 
directa del Principio de Fermat, o del recorrido de tiempo minimo. (Cantos, 
1980). 

Sea la Figura 39 en donde se han considerado dos refractores. Por el principio 
de reciprocidad la prolongacion de las dromocronas de un mismo refractor 
generadas por puntos de disparo conjugados, por ejemplo dromocrona H-l y K-J, 
deben interceptar el eje del tiempo en un mismo valor, es decir T 1 d=T la. 
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Figura 38. Ley de las velocidades aparentes. 



Figura 39. Principio de reciprocidad 
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8.6.3 Principio del intercepto en el origen. 


Sea un refractor inclinado (Figura 40) y supongase que desde la superficie se 
efectuan dos tiros en O y que se registrar! las ondas en dos puntos A y D ; el 
principio del tiempo de intercepto en el origen dice que si se prolongan las 
dromocronas a y d hasta que cortan el eje de tiempo Ot los tiempos de 
intercepto en el origen son iguales (Cantos, 1973). 



Figura 40. Tiempo de intercepto en el origen. 

8.6.4 Principio de paralelismo. 

Este principio permite deducir las dromocronas relacionadas con puntos de 
disparo intermedios entre dos puntos de disparo en los extremos del tendido; o 
deducir una dromocrona con un punto de disparo por fuera (alejado) de dos 
puntos de disparo complementarios (normal y reverso). Este principio funciona 
idealmente para refractores horizontales o inclinados pianos; en el caso de 
refractores concavos, por ejemplo hacia arriba para distancias cortas las ondas 


115 




viajaran por la superficie del refractor, pero para distancias largas viajaran (mas 
rapido) por dentro del refractor. (Siguiendo el Principio de Fermat). 

Segun la Figura 41a) la dromocrona intermedia O x HE se puede deducir por 
simple paralelas de las dos dromocronas complementarias ODE y O'FG a partir 
del punto de supuesto disparo 0 ] (Cantos, 1973). 



Figura 41. Aplicacion del principio de paralelismo. 
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8.7 METODOLOGIAS DE CALCULO 

De las suposiciones generales del metodo y de los principios y leyes se derivan 
los metodos de analisis de las dromocronas identificadas. Los metodos mas 
comunes son: tiempos de intercepto, velocidades aparentes, frentes de onda, 
tiempos de retardo y trazado de rayos. 
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Estos metodos estan completamente desarrollados y documentados en diversos 
libros de geofisica aplicada (p. ej. Jakosky, 1950; Cantos, 1973; Masuda, 1983; 
Sarria, 1996), y algunos de ellos han sido implementados en diversos programas 
de computador para agilizar los calculos. 

En general, los metodos de interpretacion y calculo en refraccion sismica se 
diferencian por las suposiciones que implican y los modelos de suelo que 
estudian, asi como por sus capacidades de aproximarse mejor a la topografia 
real de los refractores y al perfil de velocidades. 

A continuacion se describen brevemente los metodos. 

8.7.1 Tiempos de intercepto. 

Las ondas originadas en el punto de disparo una vez refractadas en los 
contactos de los medios, determinan los tiempos de llegada de las ondas en los 
geofonos en superficie. Estos tiempos se incrementan con la distancia y la 
profundidad de penetracion de las ondas. De la curva t-x, el metodo utiliza la 
pendiente de las dromocronas para calcular la velocidad de los refractores; y los 
tiempos de intercepto de las dromocronas con el eje del tiempo para calcular las 
profundidades. Este metodo es usado para modelos de un refractor piano o 
multiples refractores pianos. 

8.7.2 Velocidades aparentes. 

Basandose en el principio de las velocidades aparentes, este metodo permite la 
identificacion de velocidades y profundidades para modelos con capas inclinadas 
paralelas, utiliza los tiempos de intercepto en el origen de un tendido directo y de 
su reverso. 

El metodo supone que la velocidad de cada uno de los estratos es constante 
(medios homogeneos) y que la pendiente del refractor tambien lo es. 
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8.7.3 Frentes de onda. 


El metodo de los frentes de onda es una extension del Principio de Huygens. 
Este metodo es muy apropiado para describir refractores ondulados. 

Es un metodo grafico que se apoya en el trazado de frentes de onda 
provenientes de puntos de disparo conjugados de tal manera que se pueda 
definir un punto intermedio en el refractor, tal que la suma de los tiempos de 
viaje, entre los puntos de disparo hasta los puntos de emergencia de las ondas 
refractadas desde el punto intermedio, sea igual al tiempo total de viaje entre los 
puntos de disparo conjugados. 

8.7.4 Tiempos de retardo. 

Este es un metodo util en el caso de refractores con topografias que no son muy 
accidentadas, o que presentan muchas curvaturas. 

Utiliza la definicion del tiempo de retardo como la diferencia entre el tiempo que 
requiere la onda para recorrer la trayectoria ente el punto de disparo y el 
refractor, con el angulo critico dentro del medio superior y con su propia 
velocidad, vl, y el tiempo requerido por las misma onda para recorrer la 
proyeccion de esa misma trayectoria, pero con velocidad de refractor v2. 

8.7.5 Trazado de rayos. 

Para modelos geologicos complejos, de multiples capas con refractores de 
topografia irregular, el uso de las computadoras para el calculo y analisis de las 
curvas t-x se ha vuelto estrategico en los ultimos anos. 

El metodo mas implementado en los programas especializados de computador 
es el trazado de rayos. Este es un metodo interactive e iterativo que debe partir 
de un modelo de capas - velocidades y profundidades determinadas - calculado 
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a partir de algunos de los metodos mencionados anteriormente - por ejemplo 
tiempo de retardo para empezar a calcular diferentes modelos a partir del 
trazado de rayos sismicos desde los puntos de disparo. 

8.7.6 Cuadro comparative. 

En la siguiente pagina (Tabla 11), se incluye una comparacion de los metodos 
anteriormente descritos. 

8.8 METODOS APLICADOS EN EL PROYECTO 

En este proyecto de grado, se utilizaron los metodos de Velocidades aparentes y 
Tiempos de intercepto vs “Ray front” o “Ray tracing”, debido a que se considero 
necesario disponer solamente de las profundidades promedias al refractor y de 
sus velocidades promedio. 

En § 8.8.1 y 8.8.2, a continuacion, se explican los planteamientos de cada uno 
de los metodos. En el § 8.8.3 se explican brevemente una metodologia para 
calcular las profundidades minimas de exploracion cuando no se ha encontrado 
refractor. 

8.8.1 Tiempos de intercepto. 

Las ondas originadas en el punto de disparo que se refractan en los contactos de 
los medios, o capas de suelo, determinan los tiempos de las primeras llegadas 
en los geofonos en superficie. Estos tiempos se incrementan con la distancia y la 
profundidad de penetracion; dicho incremento se emplea para determinar la 
profundidad de los estratos. 
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Tabla 11. Comparacion de metodos de calculo e interpretacion de refraccion sismica. 


Metodo d) 

Caracteristica 

Tiempos de 

intercepto 

Velocidades aparentes 

Tiempos de 

retardo 

Frentes de onda 

Principio o 

Ley 

Tiempos de intercepto 

Velocidades aparentes; 

Velocidades crecientes con la 

profundidad. 

Principio de Fermat 

Huygens; 

Velocidades aparentes. 

Suposiciones 

Medios homogeneos; 

Velocidades crecientes con la 

profundidad. 

Medios homogeneos; 

Velocidades crecientes con la 

profundidad. 

Velocidades crecientes con 

la profundidad. 

Velocidades crecientes 

con la profundidad. 

Utilization 

Capas planas 

Capas con angulos de inclinacion 

no muy diferentes 

Refractores con relieve 

poco pronunciado 

Refractores con relieve; 

refractores fallados 

Calculos 

Rapido 

Rapido para dos refractores, 

programable en hoja de calculo 

para mas. 


Lentos; requiere trazado 

de los frentes de onda en 

cada punto 

Profundidad a 

refractor (es) 

Calculable para los puntos de 

disparo 

Calculable para los puntos de 

disparo 

Calculable para todos los 

geofonos 

Calculable para todos los 

geofonos 

Velocidades 

Se calcula una velocidad para 

cada capa 

Se calcula una velocidad para 

cada capa 

Se calcula una velocidad 

en cada geofono 

Se calcula una velocidad 

en cada geofono 

Limitaciones 

Inherentes a las suposiciones 

Inherentes a las suposiciones 

Los puntos bajo la vertical 

del geofono se encuentran 

desplazados en x 

Inherentes a las 

suposiciones 


Caso de dos estratos. 

Este caso corresponde al de un medio con velocidad vl, que suprayace uno con 
velocidad v2, que corresponde a un espacio semi-infinito, separados por una 
superficie de contacto plana como se ilustra en la Figura 42. Sea i c el angulo de 
incidencia critico de las ondas que se refractan. 
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Figura 42. Tiempos de intercepto: modelo de dos estratos pianos 


Las ondas generadas en un punto de disparo al extremo del tendido, punto A 
(Figura 42) generaran unos tiempo llegada primarios segun el frente de ondas 
directas - del medio con velocidad vl- y del frente de ondas refractadas - del 
medio de velocidad v2-. La curva tiempo distancia relacionada se muestra en la 
Figura 43. El primer tramo corresponde a los tiempo de llegada de las ondas P 
directas y el segundo tramo a los tiempo de llegada de las ondas P refractadas. 

En cada tramo, la pendiente corresponde al inverso de la velocidad vl para la 
capa superior y v2 para la capa inferior. La proyeccion del primer tramo en el eje 
de las distancias se conoce como x c , que es la distancia entre el punto de 

disparo y el sitio en superficie donde emerge la primera onda refractada, y t a , es 
el tiempo transcurrido. 
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Figura 43. Tiempos de intercepto 


En el segundo tramo los tiempos de llegada podrian relacionarse con el angulo 
de incidencia critico, i c , y la velocidad de la capa en el espacio semi-infinito, v2. 


El tiempo de viaje de la onda refractada desde el punto A hasta el punto D 
(Figura 44) es el siguiente (Sarria, 1996): 

(0 a bcd = tAB + tBC + tCD 


Como la proyeccion vertical del rayo incidente critico, AB, es igual al espesor de 
la primera capa, h, y la proyeccion horizontal es igual a h*tan(i c ) y como el 
tiempo de recorrido AB seria igual al CD, entonces: 


tAB = tCD = 



vl 
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Figura 44. Metodo del intercepto: geometria capas planas. 

y el tiempo de recorrido BC seria: 

tbC _ x-2*zl*tan(Q 
v2 

donde x, es la distancia desde el punto de disparo (desde A). 


De tal manera que el recorrido total estaria dado por la siguiente formula: 

_ ^ / cos (i c ) x - 2 * zl * tan(z' c ) 

V / ABCD = ^ 1 I 

vl v2 


Aplicando la ley de Snell, se tiene que: 


"cosO' ) x 
vl v2 


Al hacer x = 0 , en la grafica de tiempo distancia (Figura 44) el tiempo seria el 
tiempo de intercepto t 0 , de tal manera que: 
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h = ^ * t 0 * vl * cos(z' c ) 

Donde Aes la profundidad por debajo del punto de disparo. 

h = * t 0 * vl * cos(sen ( V J/ 2) 

Caso de varios estratos. 


En el caso de varios estratos pianos (Figura 45), la geometria de refraccion de 
las ondas en los estratos mas profundos es un poco mas compleja, pero se basa 
en los principios expuestos para el caso de dos capas y sus resultados se 
presentan a continuacion: 
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Figura 45. Modelos para tres refractores. 
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8.8.2 Velocidades aparentes. 

Basandose en el principio de las velocidades aparentes, este metodo permite la 
identificacion de velocidades y profundidades para modelos con capas inclinadas 
paralelas, utiliza los tiempos de intercepto en el origen de un tendido directo y de 
su reverso. 


Caso de un estrato inclinado. 


Considerese un modelo de dos capas (Figura 46), la primera con velocidad vl 
que suprayace sobre una capa de extension semi-infinita. Las capas se 
encuentran separadas por un contacto inclinado piano de angulo a 2 (Cantos, 
1980). 
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Figura 46. Velocidades aparentes: modelo de un refractor. 


Para describir el contacto, o superficie de refraccion, basta con determinar las 
profundidades en dos puntos. Usualmente se calculan las profundidades debajo 
de los puntos de disparo de una linea sismica directa y una reversa. 

Las suposiciones del metodo son: que la velocidad de los estratos es constante 
(medios homogeneos) y que la pendiente del refractor tambien lo es. 


, „ AB BC 

tAD — 1 h 

vl v2 


CD 


Donde AB, BC y CD se definen asi: 


AB = 


Z 1 

cos(/ c ) 
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CD = 


Z1 Z\- AD* sens(a 2 ) 
cos (i c ) cos (i c ) 

BC = A'D'-A'B - CD'= AD * cos(a 2 ) - AB * sen{i ] ) - CD * sen(i i ) , 

como z'j = i c => sen{i 1 ) = v J / 0 entonces 


. ^ Z2? AD* cos(a 2 ) - AB*sen(L) . ... CD 

tAD = + — * send , ) + 

vl vl 1 vl 


CD 

TAD = — * (1 - sen 2 (i , ) + -*(!- sen 2 (i x )) + 


AB 

v, 


AD * cos(a 2 ) * sen(i l ) 


Z1 * cos 2 (z'i) Z\-AD*sen{a 2 ) 2 
v 1 *cos(z 1 ) vl*cos(z' 1 ) 1 


ZD * cos(a 2 ) * sewf/, ) 


_ Zi^cosQ 't) , AD 


AD 


sen(i { -a 2 ) + 2 


H * send \ -a 2 ) 

Zl*cos(z 1 ) 


La anterior es la ecuacion total del tiempo recorrido entre Ay D. Al analizarla 
para distancias x variables, que aumentan desde A hacia D, se tiene: 

Tax = —* send, -a 2 ) + 2 * Z ‘ * cos(/l) [l] 

Vl 

En el caso de un disparo en D, el tiempo de D hacia A estaria dado por la 
siguiente formula: 

Tdx = —* sen(i \ + a 2 ) + 2 * Zl * C ° S(/l) [2] 

Vl Vl 

Segun se definio anteriormente, las velocidades aparentes para el sentido 
directo - de A hacia D como para el reverso - de D hacia A - estan definidas 
por W = dx/Zr correspond iente con la tangente en el punto de analisis (A 6 D). 
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Las velocidades aparentes entonces, se pueden determinar en funcion de la V x y 
del angulo emergente en A = /, -a 2 y en D = /j+a 2 como se muestra a 
continuacion: 

V, 

Velocidad Aparente del estrato 1 en el sentido Ascendente = V Xa = 

sen(i x - a 2 ) 

vl 

Velocidad Aparente del estrato 1 en el sentido Descendente = V xd =— 

sen(i x +a 2 ) 


Ahora bien, sabiendo que los datos conocidos por interpretacion de las 
dromocronas son V x , V aX y V a2 , hay que definir el angulo critico i x y la pendiente 

del refractor a 2 en terminos de estas tres variables conocidas. Por lo tanto en las 
2 ecuaciones anteriores se tendrian solo 2 incognitas por resolver, /, y a 2 . De tal 
manera que despejando y resolviendo se tiene que: 




-i r Fi _ x V x 
sen ( ) - sen ( — 


V, 


Id 




) 


Con base en la Ley de Snell para el angulo critico, se puede definir la velocidad 
del estrato 2, en terminos de i x o de a 2 asi: 


V 2 = 


Vx 

sen(i x ) 


V *V 

V _ 2 *_^ — *cos(a,) 

V +V 

Dos datos adicionales son determinados por las dromocronas identificadas; 
estos son el tiempo de intercepto del refractor en los puntos disparo A y D, es 
decir, el tiempo ascendente t Xa y el tiempo descendente t Xd . Estos tiempos son 

iguales a los tiempos de las ecuaciones [l] y [ 2 ] (pagina anterior), cuando x = 0. 
De tal manera que despejando en cada ecuacion en terminos de Z1 y <zl se 
tiene que : 
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Z 1 = 


T *V 
f\a M 

2 * cos(z 1 ) 

T *V 

fid y 1 

2 * cos(z 1 ) 

De tal manera que para obtener las profundidades verticales al refractor desde A 
y desde D, basta con calcular las proyecciones de Z1 y zl, asi: 


H 1 = 


1 


r *v 

1 la V 1 M 

2*cos(z 1 ) cos(aj) 

T *V 1 

hi = — — — 

2*cos(z!) cos(a,) 


Caso de dos estratos inclinados. 

El planteamiento geometrico para este caso, que se ilustra en la Figura 47 es 
similar al anterior, solo que requiere mayor manejo algebraico, de tal manera que 
se presentara solamente el plantemiento general de los tiempos de llegada en 
cada punto de disparo y las soluciones. 



Figura 47. Velocidades aparentes: modelo para dos refractores inclinados. 
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El punto de partida para el calculo de dos refractores inclinados, supone haber 
calculado previamente el angulo de inclinacion de la primera capa y las 
velocidades de la primera y la segunda. Se definen entonces los tiempos de viaje 
en el sentido ascendente del refractor 2 (T la ) y en el sentido descendente del 


mismo refractor ( T 2d ), en funcion de la distancia x: 


^ x . , _ . Z, . cost a^, + Bn) + 1 Z„ , . . 

T 2 = — * sen (B ,, -a 2 ) + — * — *-=11 + — 2 * cos( i 2 ) 

V 0 V 0 cos(a 21 ) V x 


T ld = — *sen( a 21 -a 2 ) + — 


Zi * cos(a 21 + p 21 ) + 1 z 


V r 


K 


cos(p 21 ) 


+ ■ 


22 


V t 


Derivando las anteriores ecuaciones con respecto a x y aplicando la Ley de 
Snell se tiene que: 

• _ Y 21 ~*~^21 
2 2 

= Y21 -821 


sen(i 2 ) 

Al igual que en el caso de una sola capa al hacer x = 0,T 2d yT 2a ser\ar\ los 
tiempos de intercepto en el origen. De tal manera que se puede despejar Z 22 y 
z 22 y luego reemplazarlas, igualando a H2 y hi, se obtienen las profundidades 
a I refractor: 


H 2 = 


1 


cos(a 2 ) 


? * cos(a 2 i a 3 +a 2 ) 7 

'1 / \ +Z '22 
cos(a 21 ) 


h 2 


1 


cos(a 2 ) 


_ *cos(p 21 +a 3 -a 2 ) _ _ 

Z 1 / n n +Z 22 


cos(P 21 ) 
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8.8.3 Profundidad minima. 


En la literatura no existe explicitamente una metodologia para calcular la 
profundidad minima de exploracion, cuando no se ha encontrado un refractor en 
los ensayos de refraccion sismica, para esto se deben suponer: 

a) La distancia critica igual a la distancia entre el punto de disparo y el punto de 
medicion mas lejano. 

b) La velocidad de propagacion de las ondas P en el refractor supuesto. 

El desarrollo y calculo de la profundidad minima se puede hacer por el metodo 
empleado para calcular las velocidades y profundidades previamente calculadas, 
bien sea por velocidades aparentes o tiempos de intercepto. 

8.9 ESTUDIO DE LA VELOCIDAD DE LA ONDA S 

De los ensayos de refraccion sismica se pueden obtener las velocidades de 
propagacion no solamente de la onda P, sino de la onda S. Este metodo es muy 
util en areas donde existen rocas blandas, depositos aluviales o suelos blandos, 
en los cuales pueden haber grandes variaciones entre las relaciones de la 
velocidad de las ondas P y S. Obtener la velocidad de propagacion de las ondas 
S es muy importante para determinar las caracteristicas elasticas del terreno. 

El estudio o metodo de las ondas S, utiliza los mismos principios y metodos 
analiticos explicados anteriormente en este capitulo (donde se estudia el tiempo 
de llegada de las ondas P). La diferencia radica en el uso de sensores de 
componente horizontal los cuales registran mejor la llegada de las ondas SH, 
que no se registran bien en los sensores de componente vertical de analoga 
manera como los sensores de componente vertical registran mejor la llegada de 
las onda verticales de P. 
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Este metodo esta sujeto a varias restricciones y limitaciones al aplicarlo, y no 
siempre proporciona los datos requeridos, de tal forma que otros metodos 
geofisicos como el PS logging se ha difundido mayormente (p. ej. Imai, 1975). 

Las limitaciones del metodo estan relacionadas con la energia necesaria en la 
fuente para que las ondas SH se puedan registrar a grandes distancias (p. ej. 
Sarria, 1996). 

Estas limitaciones podrian ser superadas al realizar tratamiento y procesamiento 
de las senales para realzar la llegada de las ondas S con filtros de polarizacion y 
analisis de contenido frecuencial, rutinas usualmente disponibles en paquetes de 
software sismologico especializado. 
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